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There are several possible theoretically allowed non-Abelian fractional quantum Hall (FQH) states that could 
potentially be realized in one- and two- component FQH systems at total filling fraction p = n 2/3, for 
integer n. Some of these states even possess quasiparticles with non-Abelian statistics that are powerful enough 
for universal topological quantum computation, and are thus of particular interest. Here, we initiate a systematic 
numerical study, using both exact diagonalization and variational Monte Carlo, to investigate the phase diagram 
of FQH systems at total filling fraction v = n + 2/3, including in particular the possibility of the non-Abelian 
Z 4 parafermion state. \nv = 2/3 bilayers, we determine the phase diagram as a function of interlayer tunneling 
and repulsion, finding only three competing Abelian states, without the Z 4 state. On the other hand, in single¬ 
component systems at i/ = 8/3, we find that the Z 4 parafermion state has significantly higher overlap with 
the exact ground state than the Laughlin state, together with a larger gap, suggesting that the experimentally 
observed p = 8/3 state may be non-Abelian. Our results from the two complementary numerical techniques 
agree well with each other qualitatively. 

PACS numbers: 73.43.-f, 71.10.Pm 


I. INTRODUCTION 

Multi-component fractional quantum Hall (FQH) states ap¬ 
pear in a wide variety of two-dimensional electron systems 
(2DES)', such as multilayer or multisubband quantum wells , 
systems with small Zeeman energy where the electron spin 
plays an active role, and systems with multiple valley degrees 
of freedom, such as graphene^"^, silicon**, and AlAs'^’***. These 
systems offer several tunable parameters, which allow the ob¬ 
servation of rich zero temperature phase diagrams involving 
topologically distinct FQH states even at a fixed total filling 
fraction, and indeed novel FQH phases of multicomponent 
systems have been experimentally observed. Most notable 
perhaps is the observation of the so-called 331 Abelian even- 
denominator FQH state in half-filled bilayer systems 
However, in many cases, little is known about the myriad pos¬ 
sible FQH phases and phase transitions that can be experimen¬ 
tally realized in multi-component 2DES. 

Recently, motivated by the possibility of a non-Abelian 
state at p = 5/2 in GaAs quantum wells*'*’*-^, there have 
been detailed numerical studies at total filling fraction p = 
n + 1/2 (n integer) in two-component systems'®. While 
p = n + 1/2 has been studied in great detail, the problem 
atp = n + 2/3 has received comparatively less attention from 
numerical studies'^’-”*. Such systems were first studied ex¬ 
perimentally over 20 years ago, where a two-component to 
single-component phase transition was observed in monolayer 
(presumably due to spin) and bilayer systems-'’^'"*. There are 
three Abelian EQH states that can be realized atp = 2/3: the 
330 state, 

^330 = 

i<j i,3 

consisting of two decoupled 1/3 Laughlin states in each layer 


where Zi and Wi, for f = 1,..., are the complex coordi¬ 
nates of the electrons in the two layers, and here and hereafter 
we have omitted the Gaussian factor exp(— \zi\^/fP) for 
all wave functions, a pseudo-spin singlet Abelian state, here 
called the 112 state, 

^'singlet = ^LLL ]/[ \Zi - Zj\^\Wi - /T'442, ( 2 ) 

i<j 

where T^lll is the ELL projection operator, and the particle- 
hole conjugate of the 1/3 Laughlin state, referred to here as 
the 2/3 Laughlin state, 

'I'p-H ='F’lll (3) 

i<j 

where is the wave function for the p = —2 integer 

quantum Hall state. The pseudo-spin singlet 112 state can be 
easily understood within composite fermion theory as com¬ 
posite fermions filling the lowest spin-up and spin-down levels 
in a reversed effective magnetic field''^’-*’’^". Early numerical 
work on jz = 2/3 bilayers considered the overlap of model 
wave functions with the exact ground state of the system for 
N = 6 electrons on a torus*', finding these three phases in the 
two-component 2D system (for the monolayer spinful system, 
the 330 state is unlikely). 

Different theoretical studies have suggested five possible 
exotic non-Abelian EQH states can occur at p = 2/3, yet 
have not been numerically investigated (see Table I). These 
include; the Z 4 parafermion EQH state**”*'*, 

^'Z 4 =-^[*1*330], ( 4 ) 

where A is an anti-symmetrization over all electron coordi¬ 
nates, a Eibonacci state based on SU{ 3)2 Chern-Simons the- 
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TABLE I. Candidate Abelian and non-Abelian FQH states at total 
filling fraction v = 2/3. On the sphere, these states occur at different 
shifts <S = IA — Nii, where is the number of flux quanta. The 
Fibonacci state, as a single-component system, has a shift of 6; as a 
two-component system, it has a shift of 3 per layer. 


Possible States at = 2/3 

Type 

Shift, 5 

330 State 

Abelian 

3 

Pseudo-Spin Singlet 112 

- 

1 

Particle-Hole Conjugate of 1/3-Laughlin 

- 

0 

Z 4 Parafermion 

Non-Abelian 

3 

Bilayer Fibonacci 

- 

3 

Interlayer Pfaffian 

- 

3 

Intralayer Pfaffian 

- 

3 

Bonderson-Slingerland Hierarchy 

- 

4 


Qj.y35,36^ interlayer and intralayer Pfaffian states^^-^**, 


'I'interPf 




^' 221 , 


(5) 


and 


'I'interPf = Pf 



'^ 221 , 


( 6 ) 


is thus timely to revisit the = 2/3 bilayer phase diagram 
numerically and investigate the possibility of realizing these 
non-Abelian states. 

In this work, we catTy out a study of two-component FQH 
systems at total filling fraction v = n + 2/3. We analyze 
the relative stability of the three Abelian states and the non- 
Abelian ^4 state through exact diagonalization and variational 
Monte Carlo studies that also consider the inter/intralayer 
Pfaffian states. In the lowest Landau level (LLL), our results 
are consistent with the phase diagram proposed previously '' 
and we importantly find that the Z 4 state is not competitive 
relative to the other Abelian states. However, in the limit of 
large interlayer tunneling in the second Landau level (SLL), at 
V = 8/3, our results suggest that the state is preferable rel¬ 
ative to the possible Abelian states. This unexpected new find¬ 
ing suggests the already experimentally observed 8/3 FQH 
state may be the exotic Z\ non-Abelian state, rather than the 
usual Abelian Laughlin state. Given the existence of the 5/2 
FQH state in the SLL, thought to be the non-Abelian Moore- 
Read state, the possibility that the SLL 8/3 FQH state might 
also be a (different) non-Abelian state is plausible and con¬ 
sistent with the fact that the experimental 8/3 state typically 
is considerably weaker than the 5/2 state as manifested in the 
measured energy gaps"'"'^^. 

We consider the Hamiltonian describing two quantum Hall 
layers with A total spin-polarized electrons, separated by a 
distance d, with interlayer electron tunneling strength A; 


respectively, with Xi running over the coordinates of the N 
electrons in both layers with 

^'221 = ]^(Zi - Zjf{Wi - Wjf W{Zi - WjY, (7) 
i<j i,j 


and a Bonderson-SlinglerTand hierarchy state^®. It should be 
noted that the Z 4 parafermion and the intralayer Pfaffian states 
are defined only when the number of electrons N is divisible 
by 4. 

As a result, the N = 6 overlap study^' did not actually rule 
out the possibility of having stable non-Abelian phases even 
in the lowest Landau level. The Z 4 and Fibonacci states have 
been shown theoretically to exhibit continuous phase transi¬ 
tions from the 330 state"’’"'", suggesting these states might be 
stabilized nearby more conventional ones if appropriate mi¬ 
croscopic parameters are found and tuned experimentally. The 
goal of our work is to investigate numerically the possible ex¬ 
istence of exotic non-Abelian 2/3 (or generally, n + 2/3) FQH 
states in realistic 2DES. 

The Fibonacci FQH state contains the non-Abelian Fi¬ 
bonacci quasiparticle, whose braiding statistics is known to 
be powerful enough to be utilized for universal topological 
quantum computation (TQC)"". The Z4 parafermion FQH 
state is based on the SU ( 2)4 topological quantum field theory, 
which has been discovered to allow for universal TQC 
The Bonderson-Slingerland hierarchy state dl v = n + 2/3, 
and the interlayer Pfaffian state, can also be used for univer¬ 
sal TQC if realized on topologically non-trivial spaces with 
topological operations known as Dehn twists''"'’"*'’ that can be 
realized in a physically realistic experimental setup""’’"'^. It 


N T4 

H = Y.[Y,VinUK-r^\)+^ni.r{\r} 

i<j a 



( 8 ) 


where rf is the position of the electron in layer cr, and I 
is the magnetic length. The intralayer Coulomb interaction 
is given by Vintm{r) = while the interlayer interaction is 
given by 14inter(?') = ^ dielectric of the host 

semiconductor). The interlayer tunneling term is written as 
the total pseudo-spin operator, with A the interlayer tun- 

2 

neling strength in units of ^. 


II. PHASE DIAGRAM IN THE LOWEST LANDAU LEVEL 

We first consider the 1 / = 2 jS bilayer quantum phase di¬ 
agram in the LLL. The Hamiltonian Eq. ( 8 ) has two dimen¬ 
sionless parameters: d/l, the ratio of the inter and intralayer 
Coulomb interactions, and A, the ratio of the interlayer tun¬ 
neling to the intralayer Coulomb interaction. The relative sta¬ 
bility of the three Abelian states was studied through wave 
function overlap calculations for A = 6 electrons (3 per layer) 
on the torus^*. We revisit this for larger systems using exact 
diagonalization (ED) for up to A = 12 electrons in the spher¬ 
ical geometry. In this setup, states with different topological 
orders may appear at different shifts 5 = | A — A 4 ,, where 
A$ is the number of flux quanta. 

Eigures l(a)-(f) displays our numerical results for the over¬ 
laps of the model wave functions for the 330, singlet 112, and 
2/3-Laughlin states, with the exact Coulomb ground state at 
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FIG. 1. Panels (a)-(f) correspond to shifts of 5 = 3, 0, and 1 from 
left to right for = 8 electrons, (a) - (c) show the overlap with 
i&aso, 2/3 Laughlin, and the pseudo-spin singlet (112), respectively, 
(d) -(f) show the energy gap. (g) shows the overlap with Z 4 state for 
= 12 at 5 = 3. (h) shows the gap for A^ = 12 and shift 5 = 3. 
Panel (i) displays the resulting quantum phase diagram. 
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FIG. 2. Panels (a)-(f) correspond to shifts of 5 = 3, 0, and 1 from 
left to right for A'^ = 6 electrons in the lowest Landau level, (a) - 
(c) show the overlap with ^' 330 , 2/3 Laughlin, and the pseudo-spin 
singlet ( 112 ), respectively, (d) -(f) show the energy gap. 
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shifts 5 = 3, 1, and 0, respectively, together with the energy 
gaps for N — 8 . The energy gap is taken as the difference 
between the angular momentum L = 0 ground state and the 
first excited state (if the ground state has L 7 ^ 0 , the gap is 
set to zero). We can combine the energy gaps at different 
shifts into a single function 6 {d/l, A) by choosing the max¬ 
imal gap among the different shifts. We note this gap is not 
necessarily the transport gap measured experimentally but the 
gap connected to the robustness of the phase-in many cases 
they are known to be qualitatively similar. Similar results are 
obtained for N = 6 and 10 (see Figs. 2 and 3). We do not 
compute the overlap with the Z 4 parafermion state for 5 = 3 
and N = 8 for three reasons. One is the Z 4 state is a single¬ 
component state and for N — 8 electrons there is only one 
possible L — 0 state. The second reason is the Z 4 state has 
four-electron clustering properties that cause it to vanish ex¬ 
actly unless N mod 4 = 0. Hence, one must consider at 
least N = 12 electrons, see Figs. 1(g) and 1(h). The third 
reason is the gap at 5 = 3 in the single-component limit is 
significantly below the gap at 5 = 0 corresponding to the 2/3 
Laughlin state. 

We produce a phase diagram for the bilayer system taking 
into account both the overlaps and the energy gaps at differ¬ 
ent shifts. The topological order is identified from the model 


FIG. 3. Panels (a)-(f) correspond to shifts of 5 = 3, 0, and 1 from 
left to right for A = 10 electrons in the lowest Landau level, (a) - 
(c) show the overlap with 'I'aao, 2/3 Laughlin, and the pseudo-spin 
singlet (112), respectively, (d) -(f) show the energy gap. The white 
space beyond d/l > 5.5 for 5 = 0 was not calculated and therefore 
left blank. 


wavefunction with the highest overlap with the ground state, 
and its stability is characterized by the energy gap. Figure l(i) 
shows a plot of the gap function 5{d/l, A) for N = 8 and con¬ 
tour lines showing the wave function overlaps. Our results for 
A = 6,10 electrons (not shown) are consistent with Fig. l(i). 
We emphasize this approximate phase diagram matches re¬ 
markably well with Ref. 3rs determined by wave function 
overlap and topological degeneracy on the torus. 

To investigate the relative stability of the Z 4 parafermion 
state, we consider the two-component system for A = 12 
particles. In Figs. 1(g) and 1(h), we show the overlap of 
the exact ground state with the Z 4 parafermion state, together 
with the value of the energy gap at shift 5 = 3. While the 
Z 4 state has a maximum overlap of « 0.93 in the single¬ 
component limit, the Laughlin state has a much higher overlap 
of « 0.99. Furthermore, the system possesses a much larger 
energy gap at the 2/3 Laughlin shift relative to the shift of the 
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^4 parafermion state. 

To further assess the stability of the Z 4 state compared to 
the 2/3 Laughlin state, we can consider the single-component 
limit of Eq. ( 8 ), obtained for strong tunneling A 1 and 
small d/l ^ 1, i.e., a single quantum well. The smaller 
dimension of the Fock space in this limit allows us to con¬ 
sider N = 16 electrons. Here we study a particularly realistic 
model that includes Landau level (LL) mixing (parameterized 
by the ratio of the cyclotron energy to the Coulomb energy 
K = [fujjc)//d)) and finite width of the single quantum 
well (parametrized by well width Specifically, our re¬ 

alistic Hamiltonian is 

.^realistic — Pij{m) 

m idj 

Pzjkim) (9) 

m i<j<k 

where Pijim) and Pijk{m) are operators that project onto 
pairs (i, j) or triplets {i,j, k) of electrons with relative angular 
momentum m. The Hamiltonian is parameterized by two- and 
three-body pseudopoentials Vm {w/ 1 ,k) and Vm {w/I, k), 
respectively. The two-body pseudopotentials are renormal¬ 
ized by Landau level mixing corrections which include the 
effects of virtual transitions of electrons (holes) to unoccu¬ 
pied (occupied) Landau and subband levels to lowest order 
in the Landau level mixing parameter k. Landau level mix¬ 
ing also produces an emergent three-body term that explicitly 
breaks particle-hole symmetry. This Hamiltonian is described 
in great detail in Ref. 60 and was recently implemented in a 
numerical study of the FQHE at iz = 5/2'’'. In this work 
we restrict our attention to w/I < 4; wide quantum well sys¬ 
tems are often better described as bilayers. Our results for 
the LLL are displayed in Fig. 4 where the Laughlin state is 
clearly shown to be preferable. The Z 4 overlap at 5 = 3 is 
large (« 0.82) and essentially decreases monotonically with 
K and is robust to width w/l. At S — 0 the Laughlin state has 
an overlap of nearly unity (« 0.99) and is robust to k and w/l. 
Both 5 = 0 and 3 have non-zero gaps, but the gap at 5 = 0 is 
nearly three times larger than 5 = 3. Both overlaps and gaps 
are robust to varying k and w/l. Based on these results, we do 
not expect the Z 4 state in the bilayer system at = 2/3 in the 
LLL. Our conclusions for the LLL, based on ED, are further 
corroborated with variational Monte Carlo (see below). 


III. SECOND LANDAU LEVEL 


Now that we have confidently ruled out the Z 4 state in the 
LLL, we turn our attention to the second Landau level (SLL). 
Repeating the overlap and gap calculations with the SLL pseu¬ 
dopotentials, we obtain results that are quite different from 
the LLL. In particular, we find in the single-component limit 
{d/l ^ 1 and A ^ 1) the gap at shift 5 = 3 is significantly 
larger than at 5 = 0 suggesting the ground state in this limit 
might not be the 2/3 Laughlin state, but rather an alternative 
state with 5 = 3. 

Leaving the exploration of the full two-component phase 
diagram for the next section, we now take a closer look at the 
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FIG. 4. Overlap and gap calculations for A = 16 electrons in the 
LLL in the single-component limit. Top left (a) shows overlaps with 
the Z4 state, at 5 = 3; top right (b) shows overlaps with the Laughlin 
state, at 5 = 0. Lower panels, (c) and (d), show the gaps at 5 = 3 
and 5 = 0. 
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FIG. 5. Overlaps ((a) and (b)) and energy gaps ((c) and (d)) for 
A = 16 electrons in the SLL in the single-component limit. See 
caption of Fig. 4. 


SLL in the single-component limit. We use the same realis¬ 
tic model introduced earlier (but at filling 8/3), and focus on 
the competition between the Z 4 and the Laughlin states. Sur¬ 
prisingly, we find that the Z 4 state appears favored over the 
Laughlin state according to both overlap and gap calculations, 
as shown in Fig. 5. In the SLL, the overlap with the Z 4 state is 
qualitatively similar to the LLL, i.e., it is nearly 0.83—0.84 for 
small K and decreases to zero as k is increased. The Laughlin 
state at 5 = 0 has a smaller overlap of 0.64 — 0.8, increases 
with w/l, and monotonically decreases with k. Importantly, 
the gap is approximately 1.5 times larger at 5 = 3 than it is 
at 5 = 0. Our results for A = 12 electrons are qualitatively 
similar, but with quantitatively higher overlaps for both the Z 4 
and Laughlin states. Last, we note that in the limit of zero LL 
mixing (k = 0) our Hamiltonian is particle-hole symmetric 
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FIG. 6. The quantum phase diagram determined by variational ener¬ 
gies at thickness tu = 0. The contour lines depict the energy advan¬ 
tage 5E of each dominant phase. In the SLL, the Laughlin state does 
not have a clear energy advantage over the Za state. 


lows us to obtain the state and compute its variational energy, 
i.e., the expectation value of Hamiltonian for ^Z 4 , for up to 
TV = 28 electrons, well beyond the scope of exact diagonal- 
ization. We assess the relative stability of different phases 
of the two-component system by comparing their variational 
energies. To estimate the energy per particle in the thermody¬ 
namic limit, we use quadratic extrapolation in 1 /N, weighted 
by the statistical error on each data point-^. The extrapolated 
energy has an error between 10“^ and 10“^. 

The phase diagram is determined according to the wave- 
function with the lowest energy, and we characterize the phase 
stability using the energy advantage SE of the dominant wave- 
function over its closest competitor. Figure 6 shows the con¬ 
tour plots of 5E{d/l, A). In the LLL, the phase diagram de¬ 
termined by variational energies is in qualitative agreement 
with the ED results. We find that the non-Abelian Z 4 and the 
interlayer/intralayer Pfaffian states remain energetically unfa¬ 
vorable throughout the phase diagram, and the singlet 112 oc¬ 
cupies a very small corner of the parameter space. In the SLL, 
while the Halperin 330 state still dominates at large layer sep¬ 
aration d/l, the Laughlin state at large A is now much less 
stable compared with the LLL. The main competition comes 
from the non-Abelian Z 4 state. In fact, for much of the phase 
diagram, the energy difference between the two is on the same 
order as the estimated extrapolation error (< 10“^) for the Z 4 
state. This is in strong agreement with our gap and overlap 
calculations using ED, namely, that the non-Abelian Z 4 state 
is highly competitive with the Laughlin state. Incidentally, we 
also find a small region in the parameter space that favors the 
intralayer Pfaffian state, but we do not find any parameter set 
that stabilizes the interlayer Pfaffian state. 


and our results for the spin-polarized ;/ = 8/3 state should 
translate to v = 7/3 where some recent theoretical studies 
have suggested that the FQH state at 7/3 is likely in the Laugh¬ 
lin universality class^-’^^. 

IV. VARIATIONAL ENERGIES 

The discussion so far has focused on the fate of the Z 4 state 
in comparison with the three Abelian phases, without explor¬ 
ing other non-Abelian possibilities. One may also worry about 
the finite-size effect in the ED results. To address these con¬ 
cerns, we have performed variational energy calculations that 
also include the interlayer and intralayer Pfaffian states at shift 
5 = 3, for much larger system sizes. The energy expectation 
value of Eq.(8) of the three Abelian states and the Pfaffians 
are computed using Monte Carlo for up to A = 60 electrons, 
with sample size 10^. We refer to Ref. 64 for the details of the 
Monte Carlo energy calculation and Ref. 23 for the efficient 
evaluation of composite fermion wavefunctions. It turns out 
that the anti-symmetrization used to construct T'[Eq. (4)] is 
prohibitively expensive for numerical calculations, hence, we 
leverage the Jack polynomial representation of z^ to directly 
obtain its second-quantized amplitudes'^. This technique al- 


V. CONCLUSION 

Based on our exact diagonalization and Monte Carlo stud¬ 
ies, we find that = 2/3 bilayers in the LLL, in the limit 
of weak LL mixing, most likely do not realize the non- 
Abelian Z 4 parafermion state. Most remarkably, in the single¬ 
component limit of the SLL, the non-Abelian Z 4 phase may 
be favorable for the 8/3 EQH state relative to the Laughlin 
state. Indeed, previous studies of the experimentally obtained 
energy gaps of FQH states in the SLL have already indicated 
the possibility that the electron correlations are sufficiently 
different from those of the LLL and that novel exotic states 
might be realized^^’"’*'^^''’-'’’'*’*’''’^. While weak quasiparticle tun¬ 
neling experiments through a quantum point contact^*^ suggest 
that V = 8/3 is the 2/3 Laughlin state, it cannot be considered 
to be definitive yet and more experiments are necessary. The 
implication of our finding that the observed 8/3 SLL FQH 
state may be the parafermionic Z 4 non-Abelian phase is enor¬ 
mous since this state can be utilized for universal topological 
quantum computation. 

As this work was being completed, we became aware of a 
related manuscript*^ by Geraedts et al. By utilizing primarily 
the density matrix renormalization group technique, they re¬ 
ported that the interlayer Pfaffian is stabilized for a modified 
LLL interaction with a hollow core, which is very different 
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from the LLL and SLL realistic Coulomb interactions consid¬ 
ered in our work. It is an interesting open question whether the 
disparity between the present work and that of Ref. 66 is due 
to differences in techniques, differences in models, or both. 
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